Finite low-temperature entropy of some strongly frustrated quantum spin lattices 

in the vicinity of the saturation field 
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For a class of highly frustrated antiferromagnetic quantum spin lattices the ground state exhibits 
a huge degeneracy in high magnetic helds due to the existence of localized magnon states. For some 
of these spin lattices (in particular, the ID dimer-plaquette, sawtooth and kagome-like chains as 
well as the 2D kagome lattice) we calculate rigorously the ground-state entropy at the saturation 
field. We find that the ground-state entropy per site remains finite at saturation. This residual 
ground-state entropy produces a maximum in the field dependence of the isothermal entropy at 
low temperatures. By numerical calculation of the field dependence of the low-temperature entropy 
for the sawtooth chain we find that the enhancement of isothermal entropy is robust against small 
deviations in exchange constants. Moreover, the effect is most pronounced in the extreme quantum 
case of spin |. 
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Antiferromagnetically interacting quantum spin sys- 
tems on geometrically frustrated lattices have attracted 
much attention during the last yearsi^. Whereas in 
general frustration makes the eigenstates of the quantum 
spin system very complicated, it has been found recently 
that in the vicinity of the saturation field for a wide class 
of frustrated spin lattices just owing to frustration the 
ground states become quite simple. These exact ground 
states consist of independent localized magnons in a fer- 
romagnetic environment^. They lead to a macroscopic 
jump in the zero-temperature magnetisation curve just 
below saturations^ and may provide instabilities to- 
wards lattice deformations 6 . 

In the present paper we examine the low-temperature 
entropy of several highly frustrated antiferromagnetic 
spin lattices which may host independent localized 
magnons in the vicinity of the saturation field. The 
ground state of such a system at saturation exhibits a 
huge degeneracy which grows exponentially with system 
size. For some of the considered spin systems the ground- 
state degeneracy at saturation and therefore entropy can 
be calculated exactly by mapping the localized magnon 
problem onto a related lattice gas model of hard-core 
objects. The latter models have been studi ed in many 
papers over the last few decades (see Refs. F ^8l9Ul0l and 
references therein). We complete these analytical find- 
ings for the ground-state entropy by exact diagonalisa- 
tion data for the sawtooth chain of N = 8, 12, 16 sites to 
extend our conclusions to fields below the saturation and 
to nonzero temperatures. We also examine the effects of 
exchange anisotropy, different spin values s and devia- 
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FIG. 1: The dimer-plaquette chain which hosts three localized 
magnons at fat bonds (top) and the auxiliary lattice used for 
the calculation of the ground-state degeneracy at saturation 
(bottom). The localized magnons are eigenstates for large 
enough vertical bonds J3 > J|(Ji, J2)"™- 



tions from the condition on bond strengths under which 
the independent localized magnons are exact eigenstates 
for the sawtooth chain. Finally, we discuss briefly the 
possibility of experimental verification of our findings. 

We mention recent papers of Moessner and Sondhiii, 
ZhitomirskyiS. and Udagawa et ali^ having some rela- 
tion to our investigations. These authors calculate the 
ground-state degeneracy of the 2D kagome lattice carry- 
ing classical spins for certain spin configurations (up-up- 
down structureii*^ and structures obeying a "modified 
ice rule"—) by mapping the spin problem onto a dimcr- 
covering problem on the honeycomb lattice. Udagawa et 
al. use their result to explain the residual entropy of the 
kagome ice state which occurs in the spin ice compound 
Dy2Ti2 07 under a magnetic fieldi^i^. Note, however, 
that our study refers to the frustrated quantum spin lat- 
tices. 

To be specific, we consider several geometrically frus- 
trated lattices, namely, the dimer-plaquette chaini& (Fig. 
[TJ, the sawtooth chain 1718 (Fig. [2J), two kagome-like 
chain o 19 ! 20 (Figs. El El), the kagome lattice (Fig. EJ) , and 
the checkerboard (also called 2D or planar pyrochlorc) 
lattice (Fig. EJ. The ground-state and low-temperature 
properties for the Heisenberg antiferromagnet on these 
lattices are subjects of intensive discussions. We consider 







FIG. 2: The sawtooth chain which hosts three localized 
magnons at fat V parts (top) and the auxiliary lattice used 
for the calculation of the ground-state degeneracy at satu- 
ration (bottom). The localized magnons are eigenstates for 
J 2 = y/2 (1 + A) JjiA. 




FIG. 3: The kagome-like chain of Ref. 19 which hosts three lo- 
calized magnons (marked by bold diamonds) and the auxiliary 
lattice used for the calculation of the ground-state degener- 
acy at saturation. The localized magnons are eigenstates for 
exchange bonds of uniform strength—. 



FIG. 4: The kagome-like chain of Refs. IT9II20I which hosts 
two localized magnons (marked by bold hexagons) and the 
auxiliary lattice used for the calculation of the ground-state 
degeneracy at saturation. The localized magnons are eigen- 
states for J 2 = J A 



(monomers, or monomers and dimers, or hexagons, or 
squares). 

We start with the dimer-plaquette chain shown in Fig. 
^ The auxiliary lattice (Fig. ^ bottom) is a linear chain 
of TV = jN sites which may be either occupied (if a local- 
ized magnon is trapped by the corresponding fragment of 
the initial lattice) or empty (in the opposite case). Ob- 
viously, 



W = 2^ = exp[ i log2iV 
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exp (0.1732877V). (3) 



N quantum spins of length s described by the Hamilto- 
nian 
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where the first sum runs over the bonds (edges) which 
connect the sites (vertices) occupied by spins for the 
mentioned lattices, J nm > is the antiferromagnetic 
exchange constant between neighboring sites, A is the 
anisotropy parameter, h is the external magnetic field, 
and the second sum runs over all sites. We know from 
Refs. |4j|5|, that for certain values of J nm the considered 
lattices host localized magnons (see also the correspond- 
ing Figs. n _ EJ- Due to these localized magnons the 
ground state of (JIJ for the mentioned lattices at the 
saturation field hi is highly degenerate, since the en- 
ergies of n independent localized magnon states with 
11 = 1, . . . ,?i m axj fimax ~ N are exactly the same. Be- 
cause a certain local fragment of the lattice can be occu- 
pied by a magnon or not, the degeneracy of the ground 
state at saturation, W, grows exponentially with N giv- 
ing rise to a finite zero-temperature entropy per site at 
saturation 
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The counting problem associated with the ground-state 
degeneracy can be solved after mapping the lattice which 
hosts independent localized magnons onto some aux- 
iliary lattice which is occupied by hard-core objects 



Repeating these arguments for the diamond chain2i we 
arrive at the similar result, W = 2^, however, with TV = 
IN 

3 ' 1*1 

Next we consider the sawtooth chain shown in Fig. |2] 

The auxiliary chain (Fig. [3 bottom) consists of TV = 
sites which may be filled either by rigid monomers or 
by rigid dimers occupying two neighboring sites. The 
limiting behavior of W for a large lattice TV — > 00 may 
be found in Ref. Q 



W = cxp ( log TV 



exp (0.240606TV) . (4) 



The same result (0J holds for the two-leg ladder of Refs. 
I22H23I (see Fig. la of Ref. I23T) . Similarly, for the kagome- 
like chains shown in Figs. OlElwe get 



W = exp 




exp(0.160404TV) (5) 



and 



W = exp I jjr log 1 N j w exp (0.096242TV) , (6) 



correspondingly. 

Let us pass to the 2D case. Considering the kagome 
lattice (Fig. 0) we identify i) the centres of hexagons 
(which may trap magnons) as the sites of the auxiliary 
triangular lattice and ii) the hexagons carrying localized 
magnons together with the six attached triangles as the 
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FIG. 5: The kagome lattice which hosts three localized 
magnons (bold hexagons) and the auxiliary triangular lattice 
with hard hexagons used for the calculation of the ground- 
state degeneracy at saturation. The localized magnons are 
eigenstates for exchange bonds of uniform strength—. 

shaded hexagons on the triangle lattice. Now it is evident 
that the filling of the kagome lattice by localized magnons 
corresponds to the occupation of the auxiliary triangu- 
lar lattice by hard hexagons. The hard- hexagon model 
(i.e. the triangular lattice gas with nearest-neighbor ex- 
clusion) has been exactly solved 9 . In particular, for the 
number of ways of putting hard hexagons on the tri- 
angular lattice of TV — > oo sites the accurate estimate 
is£ exp (0.333242721976. . .TV). Therefore, taking into 
account the relation between the number of sites N of 
kagome lattice and the number of sites TV of the auxil- 
iary triangular lattice, N — 3TV, we get 

W wcxp (0.1110817V). (7) 

It should be noted here that the hard-hexagon model 
also arises while calculating W for the star lattic e)?' 24 , 
however, in that case N = 6TV 

Finally we consider the checkerboard lattice shown in 
Fig. The construction of an auxiliary lattice for the 




FIG. 6: The checkerboard lattice which hosts three localized 
magnons (bold squares) and the auxiliary square lattice with 
hard squares used for the estimation of the ground-state de- 
generacy at saturation. The localized magnons are eigenstates 
for exchange bonds of uniform strength—. 



calculation of the ground-state degeneracy at saturation 
is illustrated in the lower part of Fig. Each centre of 
the square which may host a localized magnon is repre- 
sented by a site of the auxiliary square lattice. Moreover, 
the square hosting a magnon together with the eight at- 
tached triangles of the checkerboard lattice is represented 
by the shaded hard square (consisting of four elemen- 
tary cells) of the auxiliary square lattice. A one-to-one 
correspondence between independent localized magnon 
configurations and shaded hard-square configuration ob- 
viously exists. As a result, we have to consider the 
square lattice gas with TV = iN sites with nearest- 
neighbor and next-nearest-neighbor exclusion. We are 
not aware of an estimate of the entropy for such a 
model. (For the hard-square model (i.e. square lattice 
gas with with only nearest-neighbor exclusion) the num- 
ber of ways of putting hard squares on the square lattice 
of TV -> oo sites equals^ exp (1.503048082475 .. .TV).) 
A simple estimate for the lower bound for W is 2"s" w 
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exp (0.086643A). 

To summarize this part, a class of frustrated quantum 
spin lattices has a huge degeneracy of the ground state at 
saturation that leads to a nonzero residual ground-state 
entropy. For some of such models the zero-temperature 
entropy at saturation h = h% can be estimated exactly. 
These values provide the "reference points" in the low- 
temperature dependence entropy S vs. field h for the 
corresponding lattices. It is remarkably that the calcula- 
tion of W is a pure combinatorial problem and therefore 
the values of the ground-state entropy at saturation are 
not sensitive to the value of anisotropy A or the value of 
spin s. 

In what follows we discuss the dependence of the en- 
tropy S on the magnetic field h for h < hi at arbi- 
trary temperatures using full exact diagonalisation of fi- 
nite spin systems. We expect that the qualitative be- 
havior is similar for all lattices considered. Here we 
focus on the sawtooth chain, because the ground-state 
degeneracy at saturation, W, is largest and the finite- 
size effects should be smallest. We have considered saw- 
tooth ch ains of N = 8, 12, 16 sites with J\ = 1, 
J 2 



a/2 (1 + A), anisotropy parameters A 

3 



1 and 



A = 0, spin lengths s = h, 1, | at several temperatures 
kT = 0.001, 0.05, 0.2, 0.5, 1. Some of our numerical 
results are shown in Figs. 

Let us discuss the obtained results. Firstly we note 
that for several magnetic fields below saturation h < h\ 
one has a two-fold or even a three-fold degeneracy of the 
energy levels leading in a finite system to a finite zero- 
temperature entropy. Correspondingly one finds in Fig. 
(upper panel) a peaked structure and moreover a plateau 
just below saturation. However, it is clearly seen in Fig. 
(upper panel) that the height of the peaks and of the 
plateau decreases with system size TV and one has 5 = 
at T = as N — *■ oo for h < h% and h > hi, only the 
peak at h = hi does not vanish. At finite temperatures 
this peak survives as a well-pronounced maximum and it 
only disappears if the temperature grows up to the order 
of the exchange constant, Fig. (lower panel). The value 
of entropy at saturation, which agrees with the analyti- 
cal prediction Q , is almost temperature independent up 
to about kT « 0.2, see Fig. Moreover, the value of 
entropy at saturation is also almost size-independent as 
it follows from data for different N. Thus, the effect of 
the independent localized magnons which yield the resid- 
ual ground-state entropy survives at finite temperatures 
kT < 0.2 producing a noticeable enhancement in the 
isothermal entropy curve at the saturation field. 

By calculating results for A = and A = 1 and for 
s = 1 , | and s — \ we have checked that the maximum in 
the entropy at saturation for low temperatures is robust 
against exchange interaction anisotropy and appears also 
for larger spin values s. However, our numerical results 
suggest that the enhancement of the entropy at satu- 
ration for finite temperatures becomes less pronounced 
with increasing s. A simple reason for that could be the 
circumstance that the degeneracy at saturation does not 
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FIG. 7: Field dependence of the isothermal entropy per site 
at low temperatures (upper panel) and higher temperatures 
(lower panel) for the sawtooth chains of different length (s — 
|, A = 1, Ji = 1, J 2 =2). 



depend on spin value s, but the total number of state 
increases with s according to s . 

Concerning the experimental confirmation of the pre- 
dicted behavior of the entropy in real compounds we 
are faced with the situation that the conditions on 
bond strengths under which the independent localized 
magnons become the exact eigenstates&i^ are certainly 
not strictly fulfilled. For example, for the isotropic 
Heisenberg sawtooth chain (0) we have imposed J2 = 
2Ji, see Fig. Therefore, it is useful to discuss the 
"stability" of our conclusions against deviation from the 
perfect condition for bond strengths. For this purpose we 
examine numerically the field dependence of entropy at 
low temperatures for the s = h isotropic sawtooth chain 
of N = 16 sites with J\ = 1 and J2 = 1.9 and J2 = 2.1 
(Fig. IHJ . Evidently, the degeneracy of the ground state 
at saturation is lifted when J2 7^ 2 that immediately 
yields zero entropy at saturation at very low tempera- 
tures (long-dashed and short-dashed curves in the upper 
panel of Fig. |HJ. However, the initially degenerate en- 
ergy levels remain close to each other, if J 2 only slightly 
deviates from the perfect value 2. Therefore with in- 
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FIG. 8: Field dependence of the isothermal entropy per site 
of the sawtooth chain at very low temperature (upper panel) 
and at higher (but still low) temperature (lower panel) as J2 
deviates from \Jl (1 + A) J\. The corresponding dependence 
for a linear chain (LC, dotted curves) is also reported for 
comparison. 



creasing temperature those levels become accessible for 
the spin system and they manifest themselves in the en- 
tropy enhancement in the vicinity of saturation at low 
but nonzero temperatures. This can be nicely seen in 
the lower panel in Fig. [5] (long-dashed and short-dashed 
peaks in the vicinity of saturation) . To demonstrate that 
this enhancement is the effect of the localized magnon 
states in the considered frustrated quantum spin lattice 
we also report the field-dependent entropy of the s = \ 
isotropic linear chain of N = 16 sites (dotted curves in 
Fig. |SJ) which remains in this field region at least two 
times smaller. 

Let us remark that the ground-state degeneracy prob- 
lem of antiferromagnetic Ising lattices in the critical mag- 
netic field (i.e. at the spin- flop transition point), which 
obviously do not contain quantum fluctuations, has been 
discussed in the literature^. Thus, the exactly solv- 
able case of the antiferromagnetic Ising chain at criti- 
cal magnetic field provides another example of the low- 
temperature entropy enhancement at saturation. From 



Ref. we know that at zero temperature S = k log 

for s = |. Repeating the transfer matrix calculations 

for s — 1 and s = | we find instead S — k log 2 and 

S = k log ( respectively, that shows that the zero- 

temperature entropy at critical field depends on the spin 
value s. On the contrary, for the frustrated quantum 
spin lattices considered in the present paper the zero- 
temperature entropy at saturation does not depend on 
s. 

Finally, we should emphasize that there are other lat- 
tices which support the independent localized magnon 
states e.g., the 2D square-kagome lattice or the 3D py- 
rochlore lattice^. In these cases a rigorous result for 
the ground-state degeneracy at saturation is not avail- 
able, but for the existing huge degeneracy at saturation 
a lower bound is given by& W > 2" max where n max ~ N 
is the maximum number of localized magnons which de- 
pends on the lattice geometry. This leads to the con- 
clusion that the discussed low-temperature peculiarity of 
the entropy in the vicinity of saturation should also be 
present. Hence, the low-temperature maximum of S at 
saturation is a generic effect for strongly frustrated quan- 
tum spin lattices which may host independent localized 
magnons. 

From the experimental point of view the discussed ef- 
fect of the independent localized magnons on the low- 
temperature field dependence of the entropy in the vicin- 
ity of saturation may be of great importance. Really, al- 
though the most spectacular effect of the independent lo- 
calized magnons is a jump in the zero-temperature mag- 
netisation curve just below the saturations*^, it is prob- 
ably difficult to observe the jump at finite temperatures. 
The above discussed maximum in the entropy vs. field 
curve due to independent localized magnons is certainly 
easier accessible for experimental observation, since the 
isothermal entropy as a function of field can be obtained 
from a specific-heat measurement (see, e.g. Refs. Il4lll5ri . 
We also mention the significance of the maximum in the 
entropy vs. field curve to an enhanced magnetocaloric 
effect^. 

To summarize, we have rigorously calculated the fi- 
nite ground-state entropy at the saturation field for some 
strongly frustrated quantum spin lattices hosting local- 
ized magnons. To discuss the physical relevance of these 
results we have examined the field dependence of en- 
tropy at low temperatures for these frustrated systems. 
We have found that the independent localized magnon 
states produce a maximum in the isothermal entropy 
versus field curve in the vicinity of the saturation field 
at low temperatures. This effect is robust against small 
deviations from the condition on bond strengths under 
which the localized magnons exist. The reported behav- 
ior can manifest itself in the high-field specific heat mea- 
surements permitting to detect experimentally the inde- 
pendent localized magnons in frustrated quantum spin 
lattices. 
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